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ABSTRACT 

Understanding the properties of the hot corona is important for studying 
the accretion disks in black hole X-ray binary systems. Using the Monte-Carlo 
technique to simulate the inverse Compton scattering between photons emitted 
from the cold disk and electrons in the hot corona, we have produced two table 
models in the XSPEC format for the spherical corona case and the disk-like 
(slab) corona case. All parameters in our table models are physical properties of 
the system and can be derived from data fitting directly. Applying the models to 
broad-band spectra of the black hole candidate XTE J20 12+381 observed with 
BeppoSAX, we estimated the size of the corona and the inner radius of the disk. 
The size of the corona in this system is several tens of gravitational radius, and 
the substantial increase of the inner disk radius during the transit from hard-state 
to soft-state is not found. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: 
binaries — X-rays: individual (XTE J2012-F381) 



1. Introduction 

A black hole X-ray binary (BHXB) system consists of a black hole (a compact object 
with mass > 3Mq) and its companion star. Material can be transferred from the companion 
to the black hole either through Roche lobe overflow or via a stellar wind (accretion process). 
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Because the accreted matter must lose its angular momentum before it can be "swallowed" 
by the black hole, an accretion disk is usually formed and the gravitational potential energy 
of the accreted matter is released through radiation. The X-ray spectrum of a black hole 
X-ray binary system usually can be well fit with a two-component model: a black-body-like 
component and a power- law- like component (Tanaka & Lewin 1995). The black-body-like 
component turns off above 20 keV and is believed to be emitted from the accretion disk. The 
power-law-like component can extend up to several hundred keV, which suggests that a high 
temperature electron cloud (corona) exists above the accretion disk and the inverse Compton 
scattering between the disk photons and electrons in the corona is the main mechanism to 
produce this component. 

In the X-ray astronomy community, the multi-color disk model (MCD) (diskbb in XSPEC, 
Mitsuda et al. 1984; Makishima et al. 1986) plus a power-law (PL) model has been employed 
traditionally to fit the energy spectra of BHXBs, then the parameters of the MCD model 
are used to infer the physical parameters of the accretion system. However, this model is 
over-simplified in the following two aspects. 

First, the assumption that the power-law component extends straight to the low-energy 
limit of the spectrum is unreasonable, because the power-law component is believed to be 
produced by inverse Compton scattering. The most natural source of the seed photons is the 
thermal radiation from the accretion disk; they arc of the same origin as the black-body-like 
component. Therefore a low-energy cutoff must be present in the power-law component, 
because the seed photon distribution has a peak energy, below which there arc only a small 
amount of seed photons. Neglecting this low-energy cutoff by applying a simple power-law 
spectrum model would under-estimate the flux in the black-body-like component. 

Second, it has been a rather standard approach by many (including some authors of this 
paper) in the field to take the black-body-like fiux derived from the spectral model fitting as 
the true fiux of the accretion disk and infer the inner disk radius from it. This ignores the 
radiative transfer process in the production of the power-law component. Since the power- law 
component is likely produced by scattering the soft photons in the original black-body-like 
component from the accretion disk, each photon in the power-law component comes at the 
expense of a lost photon in the black-body-like component, even if no absorption occius in 
the Comptonization process. Therefore certain corrections are needed to infer the original 
fiux of the black-body-like component, otherwise the real fiux would be under-estimated, as 
realized by some authors (e.g., Kubota, Makishima & Ebisawa 2001), which might be the 
reason that the inner disk radius varies significantly when a BHXB transits from its hard 
state to its soft state or vice versa, like in XTE J1550-564 (Sobczak et al. 1999a), GRO 
J1655-40 (Sobczak et al. 1999b), XTE J2012+381 (Campana et al. 2002), etc.. 
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In order to estimate the true flux of the accretion disk and then to estimate the inner 
disk radius, one needs to estabhsh an intrinsic relation between the soft (low energy) photons 
and the hard (high energy) photons, then to make a radiative transfer correction to recover 
the original flux of the accretion disk, which makes it inevitable to calculate the Comptonized 
spectrum. The Comptonized spectrum has been computed by many authors using different 
approaches. 

One of the analytical approaches was carried out by assuming an optically thick non- 
relativistic plasma and then solving the Kompaneets equation either numerically or analyt- 
ically (e.g., Sunyaev & Titarchuk 1980; Payne 1980). This approach was improved later 
by Titarchuk (1994, T94 hereafter) and Titarchuk & Lyubarskij (1995) so that it could be 
applied to both non-relativistic and relativistic plasma, and was further verified by Hua & 
Titarchuk (1995, hereafter HT95) with different analytical approximations and Monte-Carlo 
simulations. 

Another approach is to calculate radiative transfer with Monte-Carlo simulation (e.g., 
Corman 1970; Loh & Garmire 1971; Pozdnyakov, Sobol' & Sunyaev 1977; Fenimore et 
al. 1982; etc.). This method has been implemented and described in detail by Pozdnyakov, 
Sobol' & Sunyaev (1977, PSS77 hereafter), under the assumption that the probability for 
a collision between a photon and an electron is independent of the electron energy and the 
scattering angle. It was then significantly improved by Gorccki & Wilczewski (1984) by 
removing the above inaccurate assumption. The works mentioned above take the Compton 
scattering process as the main process when computing the Comptonized spectra from hot 
plasma, some other authors also take into account other important physical processes such 
as bremsstrahlung, pair production (e.g., Skibo et al. 1995) and disk- reflection (e.g., Haardt 
& Maraschi 1991). 

Recently, Poutanen & Svensson (1996, hereafter PS96), by taking into account almost 
all important physical processes in the hot plasma and solving the radiative transfer equation 
iteratively for each scattering, built a new model to describe the emergent spectra from the 
plasma, which has been successfully tested with Monte-Carlo simulations (Stern et al. 1995). 

The above works have been applied to both AGNs and X-ray binaries (e.g., Iwasawa et 
al. 2004; Wardzinski et al. 2002; Petrucci et al. 2000; Miller et al. 2004). In this paper, we 
only consider the Compton scattering process in the hot plasma (corona) and use the Monte- 
Carlo simulation to compute the emergent photon spectra. We stress on building a simple, 
self-consistent model in which the soft photons and the hard photons can be intrinsically 
linked and the radiative transfer process can be taken into account automatically. Therefore, 
the physical parameters of the corona, as well as the parameters of the accretion disk, can 
be directly derived from the spectral fitting with the model. We will discuss the difference 



-4- 



between our model and the power-law model, especially in deriving the inner disk radius 
when apphed to BHXBs. 

This paper is organized as following. In section 2, our Monte-Carlo simulation is intro- 
duced, in section 3 the results of applying our table models to the source XTE J2012-I-381 
is reported, and summary and discussion are presented in section 4. 



2. Monte-Carlo Simulations and Table Models 

2.1. Assumptions and Simulations 

In our simulations, we assume that the accretion disk is an optically thick Keplerian disk 
(Mitsuda et al. 1984), and during the accretion process, the accretion rate (m) is constant 
and the gravitational potential energy loss of the accreted material is radiated away in black- 
body radiation locally. Therefore the temperature radial profile is T(r) ~ 7--3/4^ j,^ ^j-^g 
MCD model (Mitsuda et al. 1984; Makishima et al. 1986). We also assume that the corona 
above the accretion disk may take either spherical or disk geometry, the electron density 
distribution in the corona is uniform, and the electron energy distribution in the corona 
follows a thermal form as assumed in previous works (HT95, T94, PS96); the disk reflection 
is ignored, and the corona is in stationary equilibrium. 

We sample photons for a given combination of the parameters, i.e., the temperature at 
the inner boundary of the accretion disk, the size of the corona (for a spherical corona system, 
the size is defined as the value of its radius; for a disk-like corona system, it is defined as the 
vertical height), the temperature of the electrons in the corona, and the optical depth of the 
corona (defined as r = J UeCTdl, where Ue is electron density, ctt is Thompson scattering 
cross section and the integral is along the radial direction in the spherical corona case or 
along the vertical direction in disk-like corona case). Since the radiation is taken as a black 
body locally, the number of photons released from a ring on the accretion disk is, 

dNph = a*T^27Trdr 

(X ^ ' 



where 
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a — 
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c is the speed of light, k is Boltzmann's constant and h is Planck's constant. The photon 
energies follow blackbody distribution with temperature T(r) and the initial direction of a 
photon is sampled uniformly in Att solid angles. In a spherical corona system, the corona 
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covers a portion of the disk, then a photon emitted from the disk may or may not enter the 
corona, depending upon the emitting location and the initial direction. In disk-like corona 
systems the corona covers the whole disk, so a photon emitted from the disk always enters 
the corona first. 

After a photon enters the hot corona, the probability (p) for it to interact with an 
electron is determined by its free-path distance, 

dp oc e~'^'dTi, (2) 

and 

dri = UK-NUedl, (3) 

where ti is the optical depth of the free-path distance along the moving direction, (Tk-n is 
Klein-Nishina cross section (Berestetski et al. 1972), 



aK-N{x) = 27rro- 



and 
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hu, D, are the energy and the direction of the incident photon, respectively, and v is the 
electron velocity, all in the laboratory frame; ro is the classical electron radius (ro = e^/mc^) 
and mgC^ is the electron rest energy. It is obvious that (Tk-n depends on the incident photon 
energy, electron energy and their moving directions. In order to determine the photon free- 
path distance / in equation 3, we need to calculate Ck-n, therefore we need to know the 
electron energy (Ee) along a sampled direction, which usually is a function of the photon's 
free-path distance, i.e., the dependency relation for equation 3 is like, / — > ctk-n ^ ^ I. 
To simphfy this problem, we assume that the energy distribution in the corona is uniform, 
then the dependency between Eg and I is no longer needed. However, the sampling of the 
electron velocity is not trivial, because it is related to the input photon information and the 
total Compton cross section ax-Nix), 

p{v) (X N{v){l - V ■ n/c)aK-N{x), (5) 

where, N{v) is electron distribution which is further assumed to be isotropic and follows a 
thermal form dN^ oc e'~^''/^'^ ^/EldE^, T is the temperature of the corona. 

Equation 5 is a two-dimensional probability distribution, which is rather complicated 
to obtain and has been addressed in previous works (Sobol' 1974; Gorecki & Wilczewski 
1984). In our simulation, we adopt the algorithm developed by Hua (1997) to sample the 
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electron energy E^. and its moving direction Qg. Therefore (Tk-n can be calculated (the 
photon energy and moving direction are known) and the photon's free-path distance can be 
derived using equation 3. 

After each interaction, the scattered photon direction can be sampled by using the 
differential cross section (Akhiezer & Berestetski 1965), which is a function of electron 
velocity (v), the incident photon energy and direction, and the scattered photon energy. In 
the electron's rest frame, the formula is relatively simple (Klein & Nishina 1929), 

da 1 2 f hi'\ [ hu hi^i . ^ 



where is the scattering angle, hvi and hv are respectively the incident photon energy and 
the scattered photon energy in the electron's rest frame and related to each other by 

h. = ^ . (7) 

Therefore, after sampling the electron energy and direction, we transform the incident photon 
information to the electron's rest frame, and substitute the scattered photon energy {hv) in 
equation 6 with equation 7, then the scattered photon direction can be sampled. In the end, 
the scattered photon energy can be calculated easily. 

After zero or more scattering with electrons, an incident photon will escape and the 
escaped photons will form a spectrum, which depends on the inclination angle. In our 
simulation, we ignore the disk reflection, i.e., a photon is assumed to be absorbed by the 
accretion disk when it collides with the disk. 

It is worth noting that unlike the method described by Gorecki & Wilczewski (1984), 
we trace each single photon from its emitted location until it escapes from the disk-corona 
system, then record it as a single event, i.e., each recorded photon has exactly the same 
statistical weight as others. In this manner, we sacrifice the simulation efficiency to avoid 
calculating the escape probability and the scattering probabifity distribution in each scat- 
tering. Therefore our Monte-Carlo code has been sufficiently simplified and its efficiency 
obviously depends on the configuration of the disk-corona system (especially on the optical 
depth of the corona); the total computation time is also tolerable (e.g., it takes 5 minutes to 
produce 500,000 events for a spherical corona with optical depth around unity when running 
on a 1.50 GHz Intel Xeon CPU). 
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2.2. Simulation Summary and Comparison with the Previous Works 

For a spherical corona system, Figure 1 summarizes the dependency of simulated spectra 
on the input corona parameters (corona temperature Tc, opacity r, and size Rc) as well as 
on the input disk properties (inner disk temperature Tj„ and disk inclination angle 6). For 
a disk-like corona, the emergent spectra have the similar dependency except for Rc and 9, 
because in this case, Rc is defined as the vertical height of the corona and the corona always 
covers the accretion disk, therefore Rc is solely determined by r under the uniform density 
assumption. For the dependency on 6, in a spherical corona system, the flux of the soft 
component of the spectrum at different inclination angle differs by a cosine factor due to the 
projection effect, whereas the hard component is nearly isotropic except for at very high 6 
(Figure 2. a) because r does not depend on 6. Strictly speaking, the effective optical depth is 
slightly different for different 6 if we consider the seed photon distribution along the accretion 
disk and this difference is obvious when the inclination angle of the disk is high. However, in 
a disk-like corona system, in addition to the effective area due to the projection, the effective 
optical depth of the corona also depends on the viewing angles, therefore both the soft flux 
and the hard flux are related to ^ by a cosine-like factor (Figure 2.b). 

We compare the angle- averaged Comptonized spectra calculated from our simulation 
with those from the analytic formulae for thermal Comptonization by T94, and with those 
obtained from the iterative scattering method by PS96 (Figure 3). In our simulation and 
the calculation by PS96, the incident photons are from the center of a spherical corona and 
follow the Planck's formula, whereas in the computation by T94, the incident photons were 
assumed to follow the Wien form. The key parameters of the corona, and r, are also listed 
in Figure 3. The comparison implies that, for our work and the works of T94 and PS96, the 
results are fairly consistent with each other when the corona is optically thin (tc < 1); when 
the corona is optically thick (tc > 2) the differences between all the works become visible. 

2.3. Table models 

Based on the simulated spectra, we have built two XSPEC format table models^, 
one for spherical corona system and one for disk-like corona system. The table models 
consist of the following parameters: temperature at the inner boundary of the accretion 
disk {Tin), thermal electron temperature in the corona (Tc) , radial size of the spherical 
corona {R^ in unit of Rg — GM/c^, where M is the black hole mass), optical depth of 
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the corona (r), inclination angle of the accretion disk (9), and normalization parameter 
(Knorm) wliicli is added by the XSPEC automatically and related to the disk radius by 
Knorm = {{Rin/ km) / {D / lOkpc))^ , whcrc D is the distance of the source (see Appendix 
for detail). Because for the disk-like corona is redundant, it does not appear in the 
table models for the disk-like corona. In §3, we will apply these two models to BHXB 
XTE J2012+381. 



3. Application to XTE J2012+381 

The X-ray transient XTE J2012+381 was discovered with the All Sky Monitor aboard 
Rossi X-ray Timing Explorer on May 24, 1998 (Rcmillard ct al. 1998). Its spectra observed 
with ASCA consist of a soft thermal component (with temperature around 0.8 keV) and a 
hard power- law tail (with photon index around 3) (White et al. 1998), which is considered to 
be indicative of a black hole candidate (Lewin, Van Paradijs & Van Den Heuvel 1995). Five 
observations were carried out with BeppoSAX from 1998 May 28 to 1998 July 8 (Fig 4. a). 
The spectra with the narrow field instruments LEGS (Parmar et al. 1997) and MEGS (Boella 
et al. 1997) were extracted within a circular region with a radius of 8' centered at the source 
position. The background for these two instruments were obtained using blank sky obser- 
vations. HPGSPC (Manzo et al. 1997) and PDS (Frontcra et al. 1997) data were extracted 
using SAXDAS (hpproducts V3.0.0) and XAS respectively. All our spectral analysis were 
carried out in the energy range suggested by the BeppoSAX cookbook^: 0.12-4 keV for the 
LEGS, 1.65-10.5 keV for MEGS, 8-20 keV for the HPGSPG, 15-220 keV for PDS, and a 
2% systematic error were added during the fitting to account for the uncertainties in the 
cahbration. The pubhc software package XSPEC 11.0.1 was employed in our analysis. 

As usual, MGD plus a power-law and a Gaussian line model was used first to fit the 
data, which can describe the observations very well (refer to Figure 4(b)). The best fit 
parameters are similar to the previous results reported by Gampana et al. (2002). 

Our two table models plus a broad Gaussian line were then applied to the data, and 
these two models provide reasonably good fit to all the BeppoSAX observations (Figure 4(c)- 
4(d)) and the best fit parameters are reported in Table 1. Since we mainly focus our efforts 
on the continuum, the line parameters are not listed. 

According to the fitting results, the column density nearly remains constant with a mean 
central value 1.33 x 10^^ cm~^, and Tj„ is around 0.72 keV, which agree with the previous 
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report (Campana et al. 2002). Tc decreases during the BeppoSAX observations. In the 
first two observations, the upper hmit of is even beyond the parameter range (current 
range is 5 keV to 200 keV) in our table models. In the last three observations, the spectra 
are relatively soft and can be constrained tightly. If the system has a spherical corona, 
Rc may vary from several tens to one hundred times of the gravitational radii. The optical 
depth T is very small (< 0.5) for all the five observations, and the values are with a large 
uncertainty. The normalization parameter, Knorm.) which is related to the inner disk radius, 
also varies dramatically in the five observations. However, this apparent variation is caused 
by the uncertainty of 6 (see discussion below). In order to further check the variation, we 
re-normalize Knorm to the value at zero inclination angle, and plot the inner disk radius 
{Rin = V Knorm COS 9 km) for different observations (Figure 5). For comparison, we also plot 
the results from the MCD+PL model (Campana et al. 2002) in the figure. Even though 
with large uncertainties, derived from our models are consistent with being constant, 
within the error bars (90%) and the "sharp" increase of R^n derived from MCD+PL model 
does not show up in our results. 

4. Discussion 

In this work, we developed a Monte-Carlo code to simulate the scattering process be- 
tween the seed photons from an accretion disk and the electrons in the corona above the 
disk. Based on the simulation results, we build two simple table models for the spherical 
corona system and for the disk-like corona system, which can be used in the spectral fitting 
to directly derive the corona properties (like its size, opacity, and temperature) as well as 
the disk parameters (disk temperature, inclination and inner disk radius). 

The reason that the size of the spherical corona can be determined from the X-ray 
spectral fitting is that in our models the seed photons for the inverse Compton scattering 
come from the disk, on which the temperature is a function of the distance from the central 
black hole. The relatively small size of the corona of XTE J2012-I-381 obtained in this work 
indicates that the hard photons are mainly generated within a very small region near the 
central black hole. 

The opacity of the corona can be inferred from the ratio between the fiux of the soft 
component and that of the hard component, whereas the temperature of the corona can be 
determined by the high energy turn over in the spectrum. However, for XTE J2012-I-381, the 
T we obtained is rather small, indicating only a small portion of the disk photons were up- 
scattered to the high energies. When the source was in its hard spectral state (e.g., the first 
observation), we did not see an obvious turn-over in the spectrum, which may suggest that 
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the pure thermal corona is not a good approximation in this case and should be replaced 
by the thermal-powerlaw hybrid corona, a more realistic one proposed by Coppi (1999). 
When the source was in its soft spectral state (the last two observations), we obtained a low 
corona temperature T^, which suggests that a large uncertainty may exist in the inferred 
corona parameters, because the disk emission could be mixed with the corona emission and 
any slight change of the disk parameter (e.g., inclination angle) could result in apparently 
significant change of the corona parameters. Even though these changes will also cause some 
changes in the Knorm parameters of our table model, our conclusion marks are based on the 
re-normalized Knorm (to zero degree for example) which can be hardly impacted. 

The reason that the inclination angle can in principle be inferred is that the observed 
flux from the disk depends strongly upon the inclination angle (through a cosine factor), 
whereas the hard X-ray flux resulting from Compton scattering is almost isotropic for the 
case of spherical corona. For the case of disk-hke corona, the hard X-ray flux also depends 
upon the disk inclination angle. However, the inclination angle, which is also related to the 
ratio between the soft component flux and the hard component flux, is always coupled with 
the optical depth of the corona. Thus, inferring the disk inclination angle by spectral fltting 
may have some ambiguities. 

Determining the value of the inner disk radius is important for understanding the physics 
of the accretion disk and the black hole angular momentum (Zhang et al. 1997). The 
normalization parameter inferred with our table model, which is proportional to the inner 
disk radius squared, is nearly constant in the BcppoSAX observations, and the previously 
reported sharp increase of the disk radius when the source transits from the hard state to the 
soft state did not happen. This is because the corona, though optically thin in most cases, 
scatters some of the photons emitted from the disk and makes the observed soft component 
different signiflcantly from the original disk emission. Using the simple MCD-I-PL model and 
without considering the origin of the PL like component will certainly under-estimate the 
real flux of the blackbody-like component, especially when the source is in the hard state. 
The reason that we can estimate the real flux of the accretion disk is that the number of the 
photons emitted from the disk is known in our Monte-Carlo simulations. 

In our simulation process, we did not include the reflection mechanism to produce a 
reflection component in the emergent spectrum, but simply ignored a sampled photon when 
it collides with the accretion disk. Even though the ignored photons will not be reflected 
in the output spectrum, their corresponding seed photons are still counted as part of the 
original disk emission. A disk reflected spectrum is potentially important for probing the disk 
ionization parameter, the geometry, the inclination angle of the accreting system. When the 
disk reflected component is clearly visible in a spectrum (e.g., in Cyg X-1; Salvo et al. 2001), 
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the Comptonization model developed in this work plus a reflection model can be utilized 

to account for the observed total spectrum. As long as the combined model can fit the 
observed spectrum satisfactorily, the original disk flux could still be reconstructed from the 
normalization of the table model. 

Please also be noted that in this work we did not consider the vertical structure of 

the accretion disk (e.g., Shimura & Takahara 1995) and the gravitational potential is still 
assumed to be the Newtonian potential, therefore, the derived inner disk radius in this 
work might be substantially different from the "real" value. In fact, because of the high 
temperature (~ IkeV) around the inner disk region, the Compton scattering process may 
dominate over the free-free process in this region, and the local emergent spectrum might be 
approximated as a diluted blackbody spectrum rather than a real black body spectrum. The 
spectral difference due to this effect can be corrected by the spectral hardening factor, which 
is ~ 1.7 for an accreting system at about 10% accretion rate (Shimura & Takahara 1995). 
The correction for the real potential in the vicinity of the black hole (general relativity 
effects) could be done by introducing three factors rj, Jgr, Qgr- The t] accounts for the 
difference between the apparent and the intrinsic radii of the peak temperature, and the 
fcR is due to temperature change caused by the gravitational red shift and Doppler shift, 
whereas the qgr is related to the integrated flux change caused by the gravitational focusing 
and boosting. The value of these factors, depending on the system inclination and black 
hole spin, have been discussed in several papers (Shakura & Suniaev 1973; Cunningham 
1975; Ebisawa et al. 1994) and also have been tabulated by Zhang et al. (1997), which have 
been proved to be reasonably accurate (Gierlihski et al. 1999, 2001). We should point out 
that all these corrections, if performed in this study, will systematically change the values of 
the derived inner disk radii but will hardly change the conclusion that the inner disk radius 
is nearly constant (this situation has been proved when applying the same model and the 
above corrections to the six ultra- luminous X-ray sources (Wang et al. 2004)). 

In our table model, the normalization parameter is very sensitive to the inclination 
angle, which, as mentioned before, is coupled together with the optical depth. In order to 
reduce the uncertainty of the value of the normalization parameter in our table model and to 
infer the real flux of the accretion disk, an accurate value of the inclination angle is needed 
and the optical observations may give useful help on this issue. 
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A. Normalization parameter in our table models 

Starting from a multi-color disk black body, assuming the inner disk temperature is Tj„ 
and the inner disk radius is then the total number of photons released by the disk per 
unit time will be (by integrating Plank function from it!j„ to infinity), 

ph = 87ra*Tlrl (Al) 

where, a* — ^(f-)^/o°° ^zf Thus, the observed photon flux at distance D with inchnation 
angle 6 will be, 

Fph{9)ds = ph{e)dn 

= 8cosi9)a*Ti^lj^ (A2) 
= 8cosie)a*Tl{^yds. 

If considering the double sides of the accretion disk, there should be another factor 2 in 
equation (Al) and (A2). However, in our calculation, we only consider one side and it is 
self-consistent within our table model. 

In our simulation, for each configuration of the parameters, 500,000 photons are col- 
lected. Because the black hole (BH) absorption and the collision between the photons and 
the accretion disk (we ignore the disk reflection), the seed photon number Nseed > 500,000. 
In each simulation, six spectra will be produced by collecting the output photons in different 
directions from 0° to 35°, 35° to 45°, 45° to 55°, 55° to 65°, 65° to 75°, 75° to 85° degrees, 
corresponding to the inclination angles 22.7, 40.1, 49.9, 59.8, 69.6, 79.2 degrees respectively 
(by considering the solid angle and disk projection). To make a connection between the seed 
photon number in our simulation and the real disk flux, each spectrum had been scaled down 
by a factor f scale, which is determined by, 

f scale = Nseed/ ph (A3) 

To get the spectrum observed at different inclination angles, each angle- dependent spectrum 
has been scaled down by different area factors Asg to get the average flux at that inclination 
angle, which is determined by. 
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If we take ri„ = 1 km, D — 10 kpc in equation (A2) (in order to compare with diskbb in 
XSPEC), the relation between real r^n and normahzation of the table model is, 
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Fig. 1. — The angle- averaged simulated spectrum as a function of (a) corona temperature, 
(b) corona opacity, (c) corona size, (d) inner disk temperature, and emergent spectra as a 
function of the disk inclination angle (e) . 
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Fig. 2. — The flux dependency on the inchnation angle in a spherical corona system (a) and 
in a disk-hke corona system (b). The dotted line indicates a scaled cosine function and the 
dashed line indicates a constant value. The hard Bux indicates the photons with at least 
one scattering and the soft Bux indicates the photons without any scattering. The turn-off 
in the low angle end is artificial (the counting statistic is poor). 
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Fig. 3. — Angle-averaged Comptonized spectra E^*Flux from a spherical corona. The in- 
cident photons are assumed from the center of the corona. The diamond symbols: results 
from our simulation; dotted line: results from Poutanen & Svensson (1996); dashed line: 
results from Titarchuk (1994). Tc is the electron thermal temperature in the corona and r 
is the Thomson optical depth of the corona. See text for detail. 
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Fig. 4. — Panel (a) is the ASM/RXTE light curve of the black hole candidate XTE 
J2012+381 during its outburst in 1998; the dotted lines indicate the dates of BeppoSAX 
observations we have analyzed. The first observation is in the rising phase in which 
the hard component is relatively strong whereas the other four observations are in the 
soft state. Panel (b) shows the unfolded spectrum and the residuals for the first obser- 
vation modeling with MCD+PL with inter-stellar absorption and a Gaussian line. The 
best fit parameters are: Nh = 1.36 ± 0.02 x 10^^ cm-^, T^^ = 0.75 ± 0.01 keV, Kbb = 
1097 ±50, r = 2.22 ± 0.04, norm = 0.32 ± 0.05, Egauss = 6.1 ± 0.2 keV, ayauss = 1.0 ± 0.2 
keV , norrugauss = (3.5 ± 0.4) x 10^^, and is 277 with 244 degrees of freedom. Panels (c) 
and (d) are unfolded spectra and the residuals for the first observation modeling with our 
table models for a spherical corona system and for a disk-like corona system respectively. 
The best fit parameters are presented in Table 1. 
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Fig. 5. — The inner radius of the accretion disk inferred from our table models {triangle 
symbols) for a spherical corona system (a) and for a disk-like corona system (b) and those 
derived from spectral-fitting with MCD+PL model {diamond symbols; Campana et al. 2002). 
In order to show the comparison clearly, we shift the observation date slightly in the plot. 
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Table 1: Results from spectral fits with our table models. 



Date 


Nh 




Tc 


Rc 




e 






(MJD) 


(1022cm-2 


) (keV) 


(keV) 


(Rg) 


r 


n 


(km^) 


(/dof) 


Model A: 


spherical corona 














50961 


1 07-0.01 
-'-•'J ' +0.02 


790-0.001 
u. ( ^y+0.003 


20o;L 


2g;| 


n OOV — 0.028 
U.ZZ< ^0.030 


nt- — ** 
^0+21 


1452;|8 


292(/242) 


50976 


1 00-0.03 

-'-•'J'J+0.03 


'-'•'"'+0.005 


20o;^8 




U.iii^O.016 


rjQ — 

ZO^40 


1202^217 


207(/170) 


50982 


1 ,,,,-0.02 
J^-"J"J+0.03 


7r^-o oo3 

' -'-•J+0.004 


ioo;}j] 




U.U/U^0.004 




2357^5Q 


179(/195) 


50987 


. ,,,,-0.01 
^^•'-''-'+0.03 


U. (OU^O.006 


r — ** 

^•<5+0.8 


^^^ + 10 


u.uyo_^o.oi9 


59(/ia;) 


2432;!!^ 


197(/148) 


51002 


1 01-0.01 
J--Ol_|_o.o3 


7 00 - 0.003 
u. ( >jy_|_o.oo3 


"3-^+1.0 


qq-20 


o.097;ni7 


29;:^* 


1293_|_j^32 


175(/149) 



Model B: disk-like corona 



1 oc-0.02 
^^••^•3 + 0. 01 
01 -0.03 
-'^•'J^+0.02 
1 00-0.02 
1.00_|_o.o2 
1 Q4-OOI 
■•■•"5^+0.01 

-1 09-0.02 

-'-•'^'^+0.01 



0.003 
+0.007 
0.003 
0.002 
-0.002 
+0.001 
0.002 
+0.001 
-0.004 
+0.001 



^32~ 

40 

100;^^ 

io.5;;:2^ 
i3-o;;:I 



n6-o "o7 

U.iiO^Q QQ9 

n ni C-0.003 

U.UlO_^0.017 

U.UOU_|_Q 007 

999 — 0.013 
.ZZZ_|_o 021 
n 90-1-0.061 
U.zoi_|_o 085 



126G;^^i 
2957^]^3gg 
1292;30 
2126;^^ 

1 c:c:q-266 
lOOO_)_924 



50961 
50976 
50982 
50987 
51002 



0.722 
0.708 
0.713 
0.720 
0.701 



200 
200 



18 
69 



+ 16 

— 

+4 



1 Q — 

-■■0+28 
ro-12 
^"5+16 



37 



+26 



278(/243) 
208(/171) 
182(/196) 
193(/149) 
161(/150) 



**the limit is not reachable. 

"Normalization of the table model, Knorm = ((-Rm/km)/(Z>/10kpc))^. 



